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Abstract 

Fluctuation relations (FRs) are among the few existing general re¬ 
sults in non-equilibrium systems. Their verification requires the mea¬ 
surement of the total work (or entropy production) performed on a 
system. Nevertheless in many cases only a partial measurement of the 
work is possible. Here we consider FRs in dual-trap optical tweezers 
where two different forces (one per trap) are measured. With this 
setup we perform pulling experiments on single molecules by moving 
one trap relative to the other. We demonstrate that work should be 
measured using the force exerted by the trap that is moved. The force 
that is measured in the trap at rest fails to provide the full dissipation 
in the system leading to a (incorrect) work definition that does not 
satisfy the FR. The implications to single-molecule experiments and 
free energy measurements are discussed. In the case of symmetric se¬ 
tups a new work dehnition, based on differential force measurements, 
is introduced. This dehnition is best suited to measure free energies 
as it shows faster convergence of estimators. We discuss measure¬ 
ments using the (incorrect) work dehnition as an example of partial 
work measurement. We show how to infer the full work distribution 
from the partial one via the FR. The inference process does also yield 
quantitative information, e.g. the hydrodynamic drag on the dumb¬ 
bell. Results are also obtained for asymmetric dual-trap setups. We 
suggest that this kind of inference could represent a new and general 
application of FRs to extract information about irreversible processes 
in small systems. 


Significance Statement 

Fluctuation Relations (FRs) provide general results about the full work (or 
entropy production) distributions in non-equilibrium systems. However, in 
many cases the full work is not measurable and only partial work measure¬ 
ments are possible. The latter do not fulfill a FR and cannot be used to 
extract free energy differences from irreversible work measurements. We 
propose a new application of FRs to infer the full work distribution from 
partial work measurements. We prove this new type of inference using dual¬ 
trap optical tweezers where two forces (one per trap) are measured, allowing 
us to derive full and partial work distributions. We derive a set of results of 
direct interest to single molecule scientists and, more in general, to physicists 
and biophysicists. 
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Introduction 


Fluctuation Relations (FRs) are mathematical equations connecting non 
equilibrium work measurements to equilibium free energy differences. FRs, 
such as the Jarzynski Equality (JE) or the Crooks Fluctuation Relation 
(CFR) have become a valuable tool in single-molecule biophysics where they 
are used to measure folding free energies from irreversible pulling experi¬ 
ments BEI. Such measurements have been carried out with laser optical 
tweezers on different nucleic acid structures such as hairpins [31 a El i, 
G-quadruplexes [Hi and proteins 0 EHl EH [12], and with Atomic Force 
Microscopes on proteins na and bi-molecular complexes na. An impor¬ 
tant issue regarding FRs is the correct dehnition of work, which rests on the 
correct identihcation of conhgurational variables and control parameters. In 
the single-trap optical tweezers conhguration this issue has been thoroughly 
discussed dZlEHllin]. 

The situation about how to correctly measure work in small systems be¬ 
comes subtle when there are different forces applied to the system. In this 
case theory gives the prescription to correctly dehne the work (IFr) for a 
given trajectory (F): integrate the generalized force (/a) (conjugated to the 
control parameter. A) over A along F, Wp = /p f\dX. However, in some 
cases one cannot measure the proper generalized force or has limited exper¬ 
imental access to partial sources of the entropy production leading to what 
we call incorrect or partial work measurements. A remarkable example of 
this situation are dual-trap setups, mostly used as the high-resolution tool 
for single molecule studies (Fig. [^). In this case a dumbbell formed by a 
molecule tethered between two optically trapped beads is manipulated by 
moving one trap relative to the other. In this setup two different forces (one 
per trap) can be measured and at least two different work dehnitions are 
possible. In equilibrium conditions, i.e. when the traps are not moved both 
forces are equivalent: the forces acting on each bead have equal magnitude 
and opposite sign. On the contrary, in pulling experiments, where one trap 
is at rest (with respect to water) while the other is moved, the two forces 
become inequivalent. This is so because the center of mass of the dumbbell 
drifts and the beads are affected by different viscous drags (purple arrows in 
Fig. |If3), In such conditions theory prescribes that the full thermodynamic 
work must be dehned on the force measured at the moving trap whereas the 
force measured in the trap at rest (with respect to water) leads to a partial 
work measurement which, as we show below, entails a systematic error in free 
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energy estimates. The difference between both works equals the dissipation 
by the center of mass of the dumbbell, which is only correctly accounted for 
in the correct work dehnition. 

In this paper we combine theory and experiments in a dual-trap setup to 
demonstrate several results. First, we show that if the wrong work dehnition 
is used, free energy estimates will be hawed. The error is especially severe 
in the case of unidirectional work estimates (e.g. JE) while it inhuences bi¬ 
directional estimates to a lesser extent. This fact is not purely academical: 
measuring the force in the trap at rest is experimentally easier in dual-trap 
setups and, in fact, many groups choose to do so [31201 ED- For example 
if diherent lasers are used for trapping and detection, measuring the force 
in the moving trap poses the additional challenge of keeping the trapping 
and detection lasers aligned while moving them. Second, we show how it 
is possible, by using the CFR, to infer the the full work distribution from 
partial work measurements. We demonstrate this new type of inference in 
our dual-trap setup by showing how to reconstruct the correct work dis¬ 
tribution (i.e the one we would have measured in the moving optical trap) 
from partial work measurements in the wrong optical trap (i.e. the one at 
rest with respect to water) by using the CFR. In particular, for symmet¬ 
ric setups the correct work distribution can be directly inferred by simply 
shifting the partial work distribution. In asymmetric setups inference is still 
possible in the framework of a Gaussian approximation, but the knowledge 
of some equilibrium properties of the system is still required. This type of 
inference should be seen an example of a more general application of FRs 
which aims at extracting information about the total entropy production of a 
nonequilibrium system from partial entropy production measurements. This 
allows us to determine the average power dissipated by the center of mass of 
the dumbbell, from which we can extract the corresponding hydrodynamic 
coefficient thereby avoiding direct hydrodynamic measurements. Moreover, 
we argue this type of inference might find applicability to future biophys¬ 
ical experiments where the sources of entropy production are not directly 
measurable, e.g. ATP-dependent motor translocation where the hydrolysis 
reaction cycle cannot be followed one ATP at a time. Finally we show how, 
in symmetric setups, the work dehnition satisfying the CFR is not unique. 
In particular the differential work, based on differential force measurements 
[22], still satishes the CFR and leads to the least biased free energy estimates. 
The distinguishing feature of this work dehnition is that it completely hlters 
out the dissipation due to the motion of the center of mass of the dumbbell. 
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Throughout the paper, and for simplicity and pedagogical reasons, most of 
the derivations, calculations and experiments are shown for symmetric dual¬ 
trap setups whereas the asymmetric case is discussed towards the end of the 
paper. 


The model 

In dual-trap setups a molecule is stretched by two optical traps, the control 
parameter A being the trap-to-trap distance. Let A,B denote the two optical 
traps. When one trap (say trap A) is moved with respect to the bath while 
trap B is at rest, suitable conhgurational variables are the positions of the 
beads, both measured from the center of the trap at rest (trap B). We shall 
denote these variables by yA and ys, Fig. 33 . The total energy of the system 
is composed of three terms: 

U{yA,yB,X) = Um{yA-yB) + YyB + Y^^~y^">^ ’ 

where the quadratic terms model the potential of the optical trap and Um 
describes the properties of the tether. However one could measure the posi¬ 
tions of the beads and the trap-to-trap distance in the moving frame of trap 
A {xa, xb and —A in Fig. |D3). and the potential energy in Eq. ([^ would 
be written as: 

U\xA,XB,X) = Um{xA-XB) + Y^A+Y^^^^^'>'^ ’ 


where yA — xa = yB — xb = A {xa and xb are negative). Central to our anal¬ 
ysis will be the equation connecting the potential U and the work performed 
on the system by changing the control parameter [T] : 




ptf 

' dtdtU. 

u 

(3) 

From Eq. ([^ we get: 


r^f 

1 fAXdt, 

(4) 

with fA = -kA{yA - A). 

Inserting U' instead of 17 in (|^ gives: 


W 

= / dtdtU' = — fB^dt, 

V 't'i t/ ^2 

(5) 
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Despite of their similarity we will show that W and W are remarkably dif¬ 
ferent. In fact, from the reference frame of trap A, the bath is seen to flow 
with velocity —A (Fig. [^). Becanse of this flow an experiment in which 
trap A is moved is not Galilean eqnivalent to one in which trap B is moved. 
In the presence of a flow the connection between potential and work, Eq. 
(|^, is not valid anymore. In fact thermodynamic work measnrements mnst 
be based on the force measnred in the trap being moved. This fact has 
been discnssed in [23], yet the implications to single-molecnle experiments 
have never been pointed ont. Snmmarizing, if trap A is moved and B is at 
rest with respect to water, nsing the work W (Eq. (§) in the JE leads to 
correct free energy estimates whereas using W (Eq. (|^) in the JE leads 
to a systematic error. Below we quantify such error in detail. The dif¬ 
ference between W and W can be readily discussed in symmetric setups 
{kA = ks = k) where calculations are much simpler. To do this we switch 
to a new coordinate system: = \{yA + l/s), = tja — Ub- Here is 

the position of the geometric center of the dumbbell, while a;_ is the differ¬ 
ential coordinate ra- In this coordinate system the potential Eq. ([^ reads 
U (a:_, a:+; A) = U~{x-; A) -|- U~^{x+; A): 


U {x_,X) = Um{x-) + ^{x_ - Xf , 
f/+(x+. A) = k(x+-^^ . 


( 6 ) 


The potential energy term 11+ associated to x+ is that of a moving trap, a 
problem that has been addressed both with experiments and theory [211ES], 
while the potential energy term U- associated to X- corresponds to pulling 
experiments performed using a single trap and a hxed point. The dumbbell 
is in contact with an isothermal bath where the equilibrium state is described 
by the Boltzmann distribution and the corresponding partition function (we 
assume a weak system-environment coupling, a situation satished in our ex¬ 
perimental conditions. Section S2 in the SI). Consequently, the two degrees 
of freedom are uncoupled and the total partition function for the system 
factorizes: 


Z{X) = Z+{X)Z-{X) 

Z^(X) = j (Jxi exp (—, 
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with f3 = {KbT)~^, T being the temperature and Kb being Boltzmann con¬ 
stant. As a consequence free energy changes in the system can be decomposed 
into two contributions: 

AG' = AG'+ + AG-, = -(5-Hog Z^. (8) 

Work can also be decomposed into two contributions, each regarding one of 
the subsystems: = f dtl/^(x±; X)dt. Here 1V“ contains the work done in 

stretching the molecule while W~^ is pure dissipation due to the movement 
of the center of mass of the dumbbell. Note that: 


w = W- + W+ , w' = W--W+ , 


(9) 


which shows that the difference between W and W is entirely due to W^. 
The JE holds for W, the standard work dehnition, so that 


AG = -/3-^ log(exp i-^W)). 


( 10 ) 


Inserting Eqs. (|^,@ in (10) we get: AG+-|-AG“ = —/3“^(exp (—/3 (iy+-f hE“))). 
In symmetric setups and W- are independently distributed random vari¬ 
ables (see Sections SI and S3 in the SI) and we can conclude that: 


AG^ = -p-Hog{e-f^^^) 


( 11 ) 


Using the JE on both W and W we get two different free energy estimates: 
AG (Eq. (^) and AG' = —log(exp (—/JJU')). The error S committed 


by using W instead of W can be quantihed as: 


£ = AG-AG' = -(5-^ log 


(exp(-/jW)) 

(exp(-/?lU'))' 


( 12 ) 


From Eqs. 0 ,( [l0| ) and again using the fact and W are independently 
distributed random variables we get (exp(—/91U)) = (exp (—+ 1U“))) = 
(exp(—/31U+))(exp(—/31U“)) and similarly for W. As a consequence 

(exp(w+!) ® reduced to 


£ =-15-Hog 


(exp(-/?lU+)) 

{exp{+/5W+))' 


(13) 
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Since a;+ is subject to a quadratic potential (Eq. (|^), we expect IV~^ to be 
a Gaussian random variable. This is not true in general for W~ given that 
X- feels the nonlinear term l/m- For Gaussian Random Variables we have: 


r' log(exp(±/3W+)) = (1V+) ± (14) 

where by we denote the variance of W~^. Moreover dragging a trapped 
bead in a fluid causes no free energy change, so that: 


AG+ = -(3-^ log(exp(-/31V+)) = (1V+) -^al = 0 


(15) 


or (1V+) = 


2 " + - 


Inserting Eqs. (14) and (15) in Eq. (13) we get: 


£ = {W+} + Gl = 2{H’+>. 


( 16 ) 


Equation (16) gives £ = > 0 showing that AG' is lower than AG (Eq. 


(12)). Interestingly enough, using W instead of W in the JE leads to free 
energy estimates in apparent violation of the second law. The error on free 
energy estimates obtained using W instead of W is proportional to the mean 
work performed on the center of the dumbbell. This mean work (W^) is just 
the mean friction force times the total trap displacement AA: 


{1V+)=7+Aa , (17) 

where 7+ is the friction coefficient of the drag force opposing the movement of 
the geometric center of the dumbbell. The value of 7+ can be independently 
obtained from equilibrium measurements [26] (Sections S4,S5 in the SI). 


Differential Work Measurements 


Equations ([^ and (15) show that free energy estimates based on the standard 
work W and the differential work W~ are equivalent: 


AG = -13-^ log(exp(-/31E)) = 

= -/3-^ log(exp(-/31E-)) - log(exp(-/?lE+)) = 
= -/3-Gog(exp(-/31E-)) = AG-. 
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(18) 








We stress that this is only true for symmetric setups where W~^ and W~ are 
independent random variables. The case of asymmetric setups is discussed 
further below. Therefore W~ can be used for free energy determination, as 
it has been done in HD, although without discussion. Equation ( [l^ does 
only hold when the number of work measurements, N, tends to inhnity. 
In all practical cases we deal with hnite N and the Jarzynski estimator is 
biased [271 EH]. The bias is strongly linked to the typical dissipation Dtyp 
and a reliable estimate of free energy differences requires a number of work 
measurements which scales as W ~ exp (T*typ) [HHl [30|, so that even a small 
reduction in T*typ entails a considerable improvement in the convergence of 
free energy estimators. Moreover the bias is superadditive. Let us consider 


GWD 


m 


for simplicity Gaussian Work Distributions. In this case the bias, 
the large N limit is a function of the variance of the distribution and of 
N [27]: 

exp - 1) 


B 


GWD 

N 


2/3N 


(19) 


is a convex function of cr, and is superadditive i.e. + 0^) > 

-h B^^^((/)^). This means that, should the work be the sum of two 
independent Gaussian contributions, the bias on the sum is greater than the 
sum of the biases. Although Eq. (19) was derived under strong assumptions, 
superadditivity does also hold for other theoretical expressions for the bias 
and has been checked in our experimental data (see below). Let us introduce 
the following Jarzynski estimators for hnite N: 


1 ^ 

( 20 ) 

i=l 

AG±=-r'log-^e-'’"'* (21) 

i=l 

and the corresponding bias functions: 

Bn = AGn-AG ( 22 ) 

B^ = AG%-AG. (23) 

Since W = W~^ + W~ , superadditivity guarantees: 

Bn>B]^ + B+> B],. (24) 
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Because of Eq. (24) differential work measurements always improve the con¬ 
vergence of free energy estimates in dual-trap setups. This is especially im¬ 
portant in all those cases in which bidirectional methods (e.g. the CFR) 
cannot be used and one has to employ unidirectional methods. 


Pulling on ds-DNA 

The theory discussed so far has been put to test in a series of pulling 
experiments performed in a recently developed dual-trap optical tweezers 
setup which directly measures force in each trap [211 E2]- The setup can 
move the two optical traps independently and measure their relative posi¬ 
tion with sub-nanometer accuracy, giving direct access to both W and W. 
In these experiments 3 kb ds-DNA tethers (~ 1 pm in contour length) were 
stretched between 1 and 3 pN (Fig. [^) in a symmetric dual-trap setup 
{kA = ks = 0.02pN/nm) using 4 pm silica beads as force probes. The exper¬ 
iments were performed moving one of the two traps (trap A) with respect to 
the lab frame and leaving trap B at rest. All experiments were performed in 
PBS buffer (pH 7.4, IM NaCl). We chose cyclical protocols (A* : Aq = At^, 
where is the total duration of the cyclic protocol). The excursion of the 
control parameter, AA = At ^/2 —Aq, was varied between 200, 400 and 600 nm, 
while the pulling speed was varied between 1.35±0.05, 4.3±0.1 and 7.2±1 
pm/s. Given the force-distance curves the total dissipation along cycles was 
measured: 



The CFR [16] is a symmetry relation between the work distribution associ¬ 
ated to the forward (Pp) and time reversed (Pr) protocols: 

Pf{W) = Pr{-W) exp{/3 (IF - AG)). (26) 


In the case of cyclic protocols Pp = Pr = P and AG 
takes the form: 


P(P) =exp {^D)P{-D) 


0 so that the CFR 
(27) 


Such symmetry of the probability distribution for D can be directly tested 
in cases where negative dissipation events are observed. In Fig. we show 
measured work histograms (solid points, left hand side of Eq. J27|) and 
reconstructed histograms (open points, right hand side of Eq. (27)). If Eq. 
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(27) is fulfilled then the measured and reconstructed histograms match each 
other. A quantitative measure of the deviation from Eq. (27) can be obtained 
from the ratio P{D)/P{—D), as shown in Fig. 
that D fulhlls the FR whereas D' does not. 
pdfs of D~^ and D~, with: 


Experimental data shows 
In Fig. we show that the 


D- = 
D+ = 


D + D' 
2 

D-D' 
2 ’ 


(28) 

(29) 


are experimentally found to satisfy a FR as in Eq. (27). D~ is just the 
differential work, W~, Eq. ([^ evaluated on a cyclic protocol, whereas D^ is 
the dissipation due to the movement of the center of mass of the dumbbell. 
Summarizing, although in general W is the only observable we expect to fulhll 
a FR, in symmetric setups two new FRs emerge, for W~ and W^. In Fig. 

,B we compare the predictions of Eqs. ([T^, ( [I^ with experimental results 
for different pulling speeds. A, and different displacements AA. Equation (17) 
must be used to correct free energy estimates obtained in all those dual-trap 
setups which do not measure the force applied by the trap which is being 
moved (as in [5l [201 [21] ). The advantages of using W~ in free energy estimates 
are shown in Fig. [^. There we show the convergence of the Jarzynski 
estimator with sample size for the cycles in Fig. HP- Being evaluated over 
cycles, the expected free energy change is zero. The convergence of the 
estimator is faster for W~ than for W in the three cases (Eq. (|2^). The 
effect is enhanced in our experiments by the high pulling speed (in the range 
1-7 /im) and by the large bead radius (2 pm). Let us note that due to the 
hnite lifetime of molecular tethers and unavoidable drift effects, raising the 
pulling speed is a convenient strategy to improve the quality of free energy 
estimates. Similar results have been found also at low pulling speeds where, 
again, W' does not satisfy the CFR (Section S7 in the SI). 


Experiments on DNA hairpins 

Fluctuation theorems are used to extract folding free energies for nucleic 
acid secondary structures or proteins. We further tested the different work 
dehnitions by performing pulling experiments on a 20bp DNA hairpin (Fig. 
4A) at a 0.96 ± 0.02 pm/s pulling speed in the same dual-trap setup as in 
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the previous dsDNA experiments. In this case the work performed during 
the unfolding and refolding of the molecule were considered separately, as it 
is customary for free energy determination. In Fig. 4B we present forward 
and reverse work histograms for W,W' and W~^. Again W and W~ both 
fulhll the CFR but W shows higher dissipation than W~, resulting in slower 
convergence of unidirectional free energy estimators (Fig. 4C). The difference 
between unidirectional free energy estimates based on W and W~ is in this 
case ~ 1 KbT. As previously discussed for double-stranded DNA (Fig. 
2A), W does not fulhll the CFR and, as a consequence, unidirectional free 
energy estimates based on W are hawed. In our experimental conditions 
the error committed by using the wrong work dehnition is again positive 
and equal to T ~ 3 KbT. As previously discussed this leads to a negative 
average dissipated work, apparently violating the second law. It must be 
noted that the diherence in free energy estimates based on W~ and hF is a 
hnite-size ehect, whose magnitude decreases when an increasing number of 
work measurements is considered. On the contrary the error committed by 
using W' does not vanish by increasing the number of work measurements. 
It can be noted from Fig 4B that, although W does not fulhll the CFR and 
gives wrong unidirectional estimates, its forward and reverse distributions 
apparently cross at W = AG within the experimental error. Although this 
could be used for free energy determination the result should be taken with 
caution as we have no general proof that this should happen in all cases. 


Free energy inference from partial work mea¬ 
surements 

The JE and CFR are statements on the statistics of the total dissipation 
in irreversible thermodynamic transformations. The need to measure the 
total dissipation limits the range of applicability of these and other FRs. For 
example testing FRs concerning the dynamics of molecular motors would 
need the simultaneous measurement of both the work performed by the motor 
and the number of hydrolyzed ATPs. Here we demonstrate that, at least in 
some cases, a different approach is possible. Let us start by considering the 
simple case of symmetric setups as developed in the previous sections. In 
the experiments discussed so far there are two sources of dissipation that we 
were able to measure and characterize separately: the motion of the dumbbell 
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and the dissipation of the differential coordinate.We already learnt that W 
satishes a FR while W does not. Moreover we know that W = {W- + W+)/2 
and W = {W- — hF+)/2, and that, being and W~ nncorrelated random 
variables, W and W have the same variance. Imagine now to have only 
partial information on the system. For example, one conld be able to measnre 
force only in the trap at rest, as many experimental setnps do. With this 
information, and in absence of any gniding principle, no statement abont the 
total entropy prodnction is possible. We will hnd snch gniding principle if 
we assnme the FR to hold for W. Knowing that W and W have the same 
variance, we jnst have to shift the work distribntion P'{W') by IF = IF' + A 
to get a new distribntion that satishes the CFR. In practice this is done 
starting from the set of IF' valnes and tnning the valne of A (Fig. i\) 
nntil P(IF) = P'(IF - A) fnlhlls the CFR (Fig. j^B). In the case of the 
hairpin, this same shifting procednre is operated for both forward and reverse 
work distribntions. Again the valne of A is tnned (Fig. |^) nntil the CFR 
symmetry is recovered (Fig. IP). The nniqne valne of A that restores the 
validity of the CFR eqnals the average work dissipated by the motion of 
center of mass of the dnmbbell, giving the hydrodynamic coefficient 7+ via 
Eq. ( p!^ . Let us note that, once the work distribution in the correct trap (i.e. 
the moving trap) has been recovered, then we could also extract the correct 
free energy difference (the value of AG, Fig. HP) and infer the distribution 
for the differential work IF_ by deconvolution. 

The extension of this analysis to the asymmetric case is more complex but 
equally interesting (Fig. 6A). The decomposition of IF, IF' in IF^ and IF“ 
(Eq. (|^) is still possible although IF’'',IF“ are not nncorrelated variables 
anymore and neither IF^ nor IF“ satisfy a FR. In this general case only IF 
satishes a FR (but not IF', IF"*", IF“). Remarkably enough, in the framework 
of a Gaussian approximation, it is still possible to infer the correct work 
distribution P(IF) out of partial work IF' measurements. The analysis is 
presented in Section S6 of the SI. In this case it is enough to know the 
trap and molecular stiffnesses km i-e. some equilibrium properties of 

the system, for a successful inference. To reconstruct P(IF) both the mean 
and the variance of P'(IF') must be changed, which can be achieved by 
doing a convolution between the P'(IF') and a normal distribution; Pa,e = 
P' *A/'(A,S) where * denotes the convolution operator and A/'(A, S) is a 
normal distribution with mean A and standard deviation S. Starting from 
a distribution P'(IF') there are inhnitely many choices of A and S which 
yield a Pa,s(IF) satisfying the CFR. Indeed, let us suppose that the pair 
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A*,S* is such that -PA*,s*(hh) satishes the CFR. Then it is easy to check 
that ^ y,*'^+ 2 <pkbT satisfy the CFR for any 0 (Fig. 6B). In this 

situation the inference cannot rest on the CFR alone. Explicit calculations 
in the Gaussian case (section S6 in SI) show that variances (n^) and means 
((...)) of P{W) and P’iW') are related by an Asymmetry Factor (AF), 


AF{kA,km,kB) 


rj^^kmikA ks) 

{W) - {W) ^ ® kAke + 2k, n) 


(30) 


which only depends on equilibrium properties such as the stiffnesses of the 
different elements (section S6.3 in the SI). Knowing the AF allows us to select 
the unique pair A, S such that AF = S^/A with Pa,-e(W) satisfying the 
CFR. The inference procedure can be described with a very simple formula 
(Section S6 in the SI). The key idea is to proceed as previously done in the 
case of symmetric setups by just shifting the mean of P'{W') by a parameter 
6 until the CFR is satished, i.e. A = h, S = 0. From the values of AF and 6 
we can reconstruct P{W) by using the formulae, 


26 ^2 _ 25AF 

2-(6AF ’ “ 2-(6AF 


(31) 


The inference procedure for an asymmetric setup is shown in Fig. 6C,6D 
for cyclic ds-DNA pulling experiments. For non-cyclic pulls (AG ^ 0) the 
procedure can be easily generalized in the line of what has been shown for 
the case of hairpin in the symmetric setup (Fig. 5C). 


Discussion 

FRs are among the few general exact results in non-equilibrium statistical 
mechanics. Their validity has been already extensively tested in different 
systems, ranging from single molecules to single electron transistors, and in 
different conditions (steady state dynamics, irreversible transformations be¬ 
tween steady states, transient nonequilibrium states). At the present stage, 
the main widespread application of FR is free energy recovery from non¬ 
equilibrium pulling experiments in the single molecule field. What we are 
presenting here is a new application of FR for inference. All FRs are state¬ 
ments about the statistics of the total entropy production in a system plus 
the environment. If some part of the entropy production is missed or inade¬ 
quately considered FRs will in general not hold. This is why, for irreversible 
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transformations between eqnilibrinm states, we have a FR for the dissipated 
work (which is the total entropy production) but not for the dissipated heat 
(which is just the entropy production in the environment). The main tenet is 
now that the violation of FRs in a given setting provides useful information: 
it is an evidence that some contribution to the total entropy production is 
being missed. We have given rigorous examples in which the violation of FRs 
can be used to characterize the missing entropy production. Remarkably, in 
our model system, one could even replace the moving trap by a moving mi¬ 
cropipette, an object lacking any measurement capability, and still infer the 
work distribution exerted by that object on the molecular system (this ex¬ 
tremely asymmetric setup would still be described by Eq.(30), with —)■ oo 
and AF = AKsThm/{kB + 2km))- These results open the exciting prospect of 
extending and applying these ideas to steady state systems, such as molecular 
motors, to extract useful information about their mechanochemical cycle. 


Conclusions 

In order to give it a clear and dehnite meaning to free energy inference we have 
discussed irreversible transformations between equilibrium states performed 
with dual-trap optical tweezers. In these experiments a molecular tether is 
attached between two beads which are manipulated with two optical traps. 
The irreversible transformation is performed by increasing the trap-to-trap 
distance at a hnite speed. In this kind of transformations the dissipated work 
equals the total entropy production leading to our first result: in pulling 
experiments work W must be dehned on the force measured in the trap 
which is moved with respect to the thermal bath. The force measured in the 
trap at rest gives rise to a work dehnition, IF', which does not satisfy the 
FR and is unsuitable to extract free energy differences. We have called W a 
partial work measurement because it misses part of the total dissipation. This 
result is of direct interest to experimentalists: many optical tweezers setups 
are designed so that they can only measure W. We have thus imagined a 
situation in which W' is measurable while W is not and asked the question: 
can we infer the distribution of W from that of W'l If the question is 
asked in full generality, without any system-specific information, the answer 
is probably negative. Knowing only the extent of violation of the FR will 
be of little use, in general some additional system specihc information will 
be needed for a successful inference. Here we discussed free energy inference 
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in the framework of a Gaussian approximation, the extent to which such 
inference is generally possible should be the subject of future studies. Let us 
summarize our main results: 

• A symmetry of the system can be crucial for the inference. For left- 
right symmetric systems (as exemplihed in our symmetric dual-trap 
setup) the P{W) can be inferred from P’{W') just by imposing that 
the former satishes the CFR. When symmetry considerations cannot 
be used, the knowledge of some equilibrium properties of the system 
may suffice to successfully guide the inference (such as the stiffnesses 
for the asymmetric setup). 

• The inference process can be used both to recover the full dissipation 
spectrum plus additional information about the hidden entropy source. 
In our specific dual-trap example W does not account for the dissipa¬ 
tion due to the movement of the center-of-mass of the dumbbell and the 
inference procedure can be seen as a method to measure the associated 
hydrodynamic drag. 

• We stress the benefits of using symmetric dumbbells in single molecule 
manipulation. In this case an alternative work dehnition, the differen¬ 
tial work W~, fulhlls the CFR and is thus suitable for free energy mea¬ 
surements. Being W~ less influenced by dissipation than W, switching 
from W to W~ ensures faster convergence of unidirectional free energy 
estimates. For asymmetric setups W~ does not satisfy a FR anymore 
(only W does) and cannot be used to extract free energy differences. 

A deep understanding of how to correctly dehne and measure thermody¬ 
namic work in small systems (a long debated question in the past 20 years) is 
not just a hue detail for experimentalists and theorists working in the single 
molecule held, but an essential question pertaining to all areas of modern 
science interested in energy transfer processes at the nanoscale. The new 
added feature of free energy inference discovered in this paper paves the way 
to apply FRs to new problems and contexts. This remains among the most 
interesting open problems in this exciting held. 
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Methods 


Buffers and DNA substrates. All experiments were performed in PBS 
Buffer IM NaCl at 25°(7; 1 mg/ml BSA was added to passivate the sur¬ 
faces and avoid nonspecific interactions. The dsDNA tether was obtained 
ligating a Ikb segment to a biotin-labeled oligo at one end and a dig-labeled 
oligo at the other end. The DNA hairpin used in the experiments has short 
(20bp) molecular handles and was synthesized by hybridization and ligating 
three different oligos. One oligo is biotin-labeled and a second is dig-labeled. 
Details of the systhesis procedure are given in [SB] . 

Optical Tweezers Assay. Measurements were performed with a highly 
stable miniaturized laser tweezers in the dual trap mode [32]. This instru¬ 
ment directly measures forces by linear momentum conservation. In all ex¬ 
periments we used silica beads with 4 pm diameter, which give a maximum 
trapping force around 20 pN. Data is acquired at 1 kHz. 
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Caption of Figure 1 

Pulling experiments with dual-trap optical tweezers. A) Force-distance 
curves in a pulling experiment on a 20 bp hairpin with our dual-trap setup. 
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Figure 1: 

A molecular tether is attached between two trapped beads. By increasing 
the distance, A, between the traps the tether is stretched or released until 
some thermally activated reaction is triggered, e.g. the unfolding or folding 
of a DNA hairpin, and detected as a force jump (black arrows). The small 
force jump (0.2 pN) is due to the low-trap stiffness of our dual-trap setup 
(~0.02 pN/nm). Inset: scheme of the hairpin with color-coded sequence 
(A/T: yellow/green, G/C: red/blue). B) Pulling experiments in a dual-trap 
setup where trap A is moved at speed A and trap B is at rest with respect to 
water. A is the control parameter, ha and ys are the conhgurational variables 
with respect to the moving trap A while xa and xb are the conhgurational 
variables with respect to the trap at rest (trap B). C) pulling curves (red 
stretching, black releasing) for a 3kb ds-DNA tether in a dual-trap setup. 
Inset: the cyclic pulling protocol used in the experiments. 
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Caption of Figure 2 


Work Distributions. Work statistics obtained on cyclic protocols with 
a 200 nm excursion with different pulling speeds (columns). Four different 
work observables are considered. In each case the solid points are direct 
work measurements (l.h.s. of Eq. ([27|)). In order to improve statistics 
these distributions are calculated as the convolution of the distribution of 
the work performed while stretching with that performed while releasing. 
For that we took all forward and reverse work values Wp, Wp and combined 
them W = Wp + Wp in order to get a joint work distribution for the cycle. 
Open symbols are the reconstructed histogram (r.h.s. of Eq. (27)). Different 
columns refer to different pulling speeds,A , as shown on top. A) Comparison 
of the measured and reconstructed distributions according to the two defini¬ 
tions of Eq. (25) The distribution for D satisfies the CFR Eq.(27), 

i.e. the measured and reconstructed distributions superimpose. The distribu¬ 
tion for D' does not satisfy it. Horizontal error bars represent the systematic 
error in work measurements, while vertical error bars denote statistical er¬ 
rors. B) The CFR, Eq. (27) is satished within the experimental error for 
D, and D~ but not for D'. C) Comparison between the measured and 
reconstructed distributions for D~ and (Eq. 


and Eq. ([M|)) both of 
which satisfy the CFR. Horizontal error bars represent the systematic error 
in work measurements, while vertical error bars denote statistical errors. 


Caption for Figure 3 


Bias in unidirectional free energy estimators. A) Error {£) on free 


energy estimates Eq. (12) committed by using JE for W. Circles, diamonds 
and squares refer to excursions of 200 nm, 400 nm and 600 nm. Three dif¬ 
ferent pulling speeds were considered in each case. Note that in this case 
work is evaluated over closed cycles (AG = 0) and the error is defined as: 
£ = —log(e“^^ ). B) {W^) displays a bilinear dependence on pulling 
speed A (Main Figure) and AA (Inset) as expected from Eq. Contin¬ 

uous lines are linear hts to the experimental results which contains a single 
htting parameter. The shaded area in the inset corresponds to the region 
within one standard deviation from the expected value of (W^) based on 
equilibrium measurements of 7+ (see Section 4,5 in the SI). C) Experimental 
bias measurements from the cycles shown in Fig. Ep. The plots show the 
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bias (Eqs. (22),(23)) as a function of the number of work measurements, N. 
The three plots correspond to different pulling speeds (7.2 pm/s, 4.3 pm/s, 
1.35 pm/s). Interestingly ^ which guarantees faster convergence of 
free energy estimates. Moreover Bj^i is also larger than the sum B~ + B~^ 
(dashed line) i.e. the bias is superadditive (cf. Eq. (24)). The error bars 
represent the statistical error on free energy determination, not including sys¬ 
tematic calibration errors in force and distance. Continuous lines show the 
theoretical predictions from Ref. [28] for Gaussian work distributions. Note 
that these are not fits but predictions which only use the mean dissipation 
as input parameter. 


Caption for Figure 4 


Work measurements on a DNA hairpin. A) Scheme of the experimen¬ 
tal setup (beads and hairpin not to scale). The hairpin is presented with 
color-coded sequence (A/T: yellow/green, G/C: red/blue). B) Work mea¬ 
surements upon unfolding and refolding according to three different work 
dehnitions: W (upper panel), W (central panel) and W~ (lower panel). In 
this experiment pulling speed was 0.96 ± 0.02 pm/s. Open symbols show an 
estimate of the refolding work distributions reconstructed from the unfold¬ 
ing distribution via the GFR (Eq. (26)). W and W~ fulhll the fluctuation 


theorem while W does not. Horizontal error bars represent the systematic 
error on work measurements, while vertical error bars denote statistical er¬ 
rors. The contribution of trap, handles and single stranded DNA have been 
removed as detailed in [B]. G) Unidirectional estimates for the free energy 
from the unfolding work distribution. The optimal estimator based on W~ 
(red) converges to the correct value AGq = 51 A^T as measured from bi¬ 
directional estimates. The estimator based on W (green) shows a larger bias 
and overestimates the free energy by ~ 1 KbT. The estimator based on W 
(blue) converges to a wrong free energy difference (AG — £) which is ~ 3 
KbT below the correct value, against the second law. Note that the error 
committed by using W is due to finite-size effects and decreases when more 
unfolding curves are measured. In contrast the error committed by using W 
remains hnite for all sample sizes. The error bars represent the statistical 
error on free energy determination and do not include the systematic error 
due to force and distance calibrations. 
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Caption for Figure 5 

Inference of P{W) from partial work measnrements in the symmet¬ 
ric case. A) The distribution P'{W') in the case of the dsDNA tether, for 
the work measured in the wrong trap, W. The distribution does not fulhll 
the CFR. To recover the correct work distribution P{W), W is shifted by a 
constant amount A. The shifted distribution is tested for the CFR by dehn- 
ing the function; iJ(hF') = log(P'(hF' — A)/P'(—hF' — A)). The prediction 
by the CFR is H{W') = (3 iW' — AG) which can be tested by determining 
the H function. B) Evolution of H{W') as a function of W for different 
values of A. The value A* for which the slope of H{W') is equal to one 
(work being measured in KbT units) determines the correct work distribu¬ 
tion P(W) (A ~ 15, inset). C) In the case of bidirectional measurements 
both the forward and the reverse work distributions P’piW') .n( —W') are 
shifted by an amount A in opposite directions. D) Evolution of the func¬ 
tion H = log(P^(lF' — A)/Pp{—W' — A)) as a function of A. Again the 
CFR predicts H should be linear in W with slope one. E) Inference of 
the correct work distributions and AG measurement. For each value of A 
a linear fit A{A)W + P(A) to H is performed. The value A* for which 
A (A*) = 1 (A* ~ 7 in the figure), is the shift needed in order to recover the 
full work distribution P{W) from partial work measurements in the wrong 
trap. Moreover the CFR implies B{A*) = —AG (AG ~ 60KbT). 


Caption for Figure 6 

Inference of P{W) from partial work measurements in the asym¬ 
metric case. (A) Equilibrium force distributions at 2pN for the 3kb ds-DNA 
tether measured in an asymmetric dual-trap setup (/ca = 0.012 pN/nm,/cB = 
0.003 pN/nm,/cm = 0.0027 pN/nm). (B) Convolution of P'(hF') with dif¬ 
ferent Gaussian distributions corresponding to different pairs A, S.We show 
P'(hF) (bine filled circles), P'(—hF) exp(/dhF) (blue empty circles), P(hF) 
(green filled circles). Among all different Pa,s only one matches the correct 
work distribntion P(hF), i.e. only one reconstrncted distribntion is physi¬ 
cally correct (rightmost graph, with A = 7.5 KbT, = 8 (P^T)^). In 
this sitnation the inference cannot rest on the CFR alone, and additional 
information is reqnired to infer P(hF). (C) Asymmetry factor (AF) as a 
fnnction of a; = kA/kB for different valnes oi y = km/kB- The blue (red) 
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circles indicate the symmetric {AF = 0) and asymmetric {AF ~ 1) cases 
respectively. (D) The AF dehned by T? = AF x A (red line) and the CFR 
invariance A = A* + 0, + 2(1)KbT for any 0 (blue line), do select 

a narrow range of possible pairs (A,S) at the intersection between the blue 
and red lines. The intersection region is compatible with the parameters (A 
=8 iF^T, YA =7.2 (iF^T)^ ) describing the true correct work distribution 
PiyV) (black point). 
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In this document we analyze in further detail several aspects of the manuscript regarding theory and 
experiments. Concerning the theory we study in detail a dynamical model describing the experiments 
reported in the main text. This model is then used to shed light on several aspects of the discussion 
in the main text. In section SI we introduce the model in the symmetric case. In section S2 we prove 
that x+ and x- are independent quantities in our model, which is a central result in our discussion. In 
section S3 we do the same for W+ and W-. These conclusions are supported by a statistical analysis 
of the experimental data. In sections S4 and S5 we comment on the definition and meaning of the 
hydrodynamic parameters of our model. One of these parameters is 7 +, which was used in estimating 
the missing contribution to the total dissipation when using the wrong work definition W'. We will 
also recall how to measure these parameters from equilibrium experiments and present measurements 
on different DNA tethers. In section S 6 we will discuss free enerfy inference in asymmetric setups and 
in section S7 we present measurements at lower pulling speeds. These measurements show that correct 
work measurements is important even at low irreversibility conditions. 


SI A MODEL FOR THE EXPERIMENTAL SETUP 


We consider a model for the experimental dual-trap setup shown in Fig. SI. In this setup two focused 
laser beams form two optical traps, which allow us to manipulate a dumbbell formed by two trapped 
beads coupled through a molecular tether. The model is one dimensional: the dynamics takes place 
along the direction of the tether, j/-axis in Fig. SI. Should the tether not be oriented in the plane 
perpendicular to the direction of the laser beams, the analysis is more involved as detailed in Ref. [1]. In 
this section we will assume the two traps are identical and harmonic, so that the total potential energy 
reads: 


CA'ot {yA,yB) 


2^1 + 2^y^~ ~ y^^ ■ 


( 1 ) 


A more general discussion, taking into account asymmetric setups, will be differed to section S 6 . Note 
that we choose the center of trap B (left trap) as the origin of our reference frame and that A is the 
position of trap A relative to trap B. A dynamical model for our system can be obtained using Langevin 
equations subject to a total potential as defined in Eq. (1): 


y = -MV17TOT(yA, J/B, At) + 77 , 


( 2 ) 


where y = {yA,yB), V = {dy^,dyg), rj — {r]A,VB) is a vector formed by the noise terms affecting 
the two beads and jj, is the mobility tensor. We assume the noise to be 6 correlated in time, and the 
Fluctuation-Dissipation Theorem to be fulfilled: 


= KBTnij6{t - s), (3) 

i.e. that the covariance of the noise is proportional to the mobility /i. In the above expression Kb is 
the Boltzmann constant and T the absolute temperature. In the vector Equation (2) the mobility is a 
tensor or matrix. Its diagonal terms describe the friction affecting the motion of the beads, while the 
off-diagonal terms describe hydrodynamic interactions between the beads: 

iA^{yA,yB) T^BiyA,yB) \ 

\ ^BAiyA,yB) lB^{yA,yB) ) ’ 

The validity of Eq. (3) for hydrodynamic fluctuations in a pulling experiment is a standard assumption, 
and it is a necessary hypothesis in the derivation of FRs such as Jarzynski equality or the Crooks 
fluctuation relation. The analysis of this equation is greatly simplihed by the following two features, 
which can be easily matched in our experimental setup: 

• The setup is fully symmetric, meaning not only kA = kB but also ta = tb where rA,f’B are the 
bead radii. 
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• The hydrodynamic coefhcients 7 , T between the beads do not significantly change in the force range 
of a pulling experiment. In section S4 we show that this condition is met for tethered molecules 
ranging from 20 nm to 1 /xm in contour length to a 10 % accuracy. 

This guarantees that the mobility matrix is symmetric and constant: 


M = 




Under these conditions, the transformation a;+ = va+vb ^ decouples the potential into two 

separate terms, 

Utot{x-,x+-, A) = U~{x--,\) + U^{x+\ A). 

Moreover the symmetry requirements cited above have the following two important consequences: 


• It diagonalizes the mobility tensor jj, in the form: 


- f 0 

0 7 _^ 


• It diagonalizes the noise covariance, i.e. it gives rise to two statistically independent noise contri¬ 
butions, and r/_ which separately affect the new coordinates x+ and X-: 


{ri±{s)r]±{t)) = KBTj±d{t - s), (? 7 +(s)r/_(t)) = 0. (4) 


In the main text we have shown how the partition functions for x- and a;+ decouple. From previous 
equations one can easily see that the full right-left symmetry of the setup guarantees that also the 
equations of motion for x~^ and x~ decouple completely: 

7+^*“^ = -dx+U^{xt, At) + vt (5) 

1-xT = -dx-U~ {xt , At) -f • (6) 


S2 THE DIFFERENTIAL X_ AND CENTER-OF-MASS A+ COORDINATES ARE 

STATISTICALLY INDEPENDENT 

The factorization of the partition function (Eq .8 in the main text) plays an important role in the deriva¬ 
tion of our results. For symmetric systems such factorization is directly related to the statistical indepen¬ 
dence of the center of mass and differential coordinates. Here we experimentally check this factorization 
by showing how we can generally describe the state of our system using a decomposition of the measured 
forces into two statistically independent coordinates. For linear systems such decomposition is always 
possible whereas for general non-linear systems such decomposition is only possible in the symmetric 
case. Imagine that the equilibrium distribution of our system is dictated by a Hamiltonian depending 
on two coordinates: U{yA,yB), which takes the form C/+(x+) -|- f7_(x_) for a general transformation, 
x- = X-{yA,yB) and a;+ = x+{yA,yB)- In this situation, a;+ and x_ are statistically independent: 

P{x+\x-) = P{x+) oc exp (-/3[/+(x+)), (7) 

and in particular have zero covariance: 

(x+x-) = {x-){x-). ( 8 ) 

A further consequence of this decomposition of the Hamiltonian is the factorization of the partition 
function discussed in the main text. For simplicity we will be considering forces instead of distances, 
as our instrument directly measures forces. The independence of the coordinates x^ and X- directly 
translates into the independence of the corresponding forces: 

/+ = -dU+{x+) /- = - 


dU-{x-) 


(9) 
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indeed: 

(/+/-) = J dx+dx-f+f-e~^^ = J dx+f+e~l^^+ j dx-f-e~^^- = (/+)(/->• ( 10 ) 

We shall use the covariance as a measure of linear dependence and a Pearson’s test as a measure of 
dependence in general. The x^ test is described below in section S3 of this document. The transformations 
we use to uncouple the two degrees of freedom are rotations of the form: 

= cos{(p)fA + sm{(p)fB (11) 

/- = - sin(()))/A + cos((())/b (12) 

We present measurements of equilibrium force fluctuations obtained using different tethers. In Figure 

52 we describe measurements performed on Ikb dsDNA tethers stretched at 2 pN. Experiments were 
performed using PBS buffer (pH 7.5), IM NaCl at 25 Celsius degrees just as in the pulling experiments 
described in the main text. The force of 2pN was chosen to fall in the range explored in the pulling 
experiments. The first set of measurements was performed using symmetric traps kA — ks — 0.02 
pN/nm. The second set of experiments was performed with asymmetric traps, kA = 0.012 pN/nm, 
ks = 0.003 pN/nm. In both cases it is possible to find a coordinate system such that the two coordinates 
are independent and in both cases this coordinate system is obtained by a rotation of the vector {/a, /b) 
(Fig. S2A,C). The difference is that in the symmetric case the new variables correspond to the center- 
of-mass and differential coordinate as discussed in the main text (Fig S2B), while they have a different 
definition (i.e. they are not generated by a tt /4 rotation) in the asymmetric case (Fig. S2D). In Figure 

53 we report similar measurements performed on a DNA hairpin (the sequence is the same reported 
in the main text) in a symmetric setup. Equilibrium traces were acquired at different forces. From 
top to bottom we can see the hairpin is I) completely folded, 2) preferentially folded, 3) preferentially 
unfolded, 4) completely unfolded. Force fluctuations are clearly non-Gaussian, a characteristic of non¬ 
linear two-state systems. Again, for the symmetric setup, our data show that the center of mass (a;+) 
and differential coordinate (x_) are independent. We note that in the more general non-linear case and 
for asymmetric setups no independent coordinates can be defined. 


S3 THE DIEEERENTIAL WORK W- AND IT+ ARE STATISTICALLY INDEPENDENT 


An immediate consequence of the decoupling of the equations of motion is the independence of IF+ and 
W-: 


W± 


d\U±(x±, X)Xdt. 


Now, since W~^ does only depend on x~^ and W~ does only depend on x~ we can conclude that the two 
quantities are independent. Of course this results holds under specihc assumptions and is, in the end, 
just a property of the model. From the experimental point of view one can test the independence of 
and W~ by statistical analysis, studying their linear (Pearson) correlations or using a Pearson’s x^ test 
for dependence. The first method is sensitive to linear correlation between two observables, while the 
second, more stringent test, is sensitive to generic correlations. The x^ test goes as follows: to test the 
independence of two random variables {A, B) from a finite number of measurements we first define the 
empirical marginals, i.e. separate histograms of the two variables: 


N 


'y ^ d{Ai € [oj, aj+i))/At 


e 


b 

j 


N 


'^5{Bi € [hj,hj+i))/N, 


(13) 


(14) 
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here 6{x € C) is 1 if a; belongs to the set C and is 0 otherwise. We then compare the actual two-dimensional 
histogram, Cij to the product of the marginals; 


X 


2 


i^ki 

paj) 


(15) 


with Cjk = S [aj,aj+i))^{Bi € \bk,hk+i))/N. The value obtained through Eq. (15) is then 

compared to the distribution. The two variables are considered independent if the value satisfies 
5% significance test. The results of the two statistical tests is reported in Table SI and S2. These tables 
compare the result of the tests performed on the pair W,W' or on the pair 1T + ,1T_, for the forward, 
the reverse and the cyclic protocols. In each case three different pulling speeds are considered. The two 
different tables refer to two different dumbbells. In most cases the pair W, W' shows higher covariance 
than the pair W+, W-. Moreover the test detects dependence between W, W but not between W+, W- 
(in the y^ column 1 denotes a positive test for dependence and 0 a negative test). 


S4 DIRECT MESUREMENT OF THE HYDRODYNAMIC PARAMETERS 


In the main text we have already shown how, in a symmetric setup, by changing to the coordinate system 
defined by the center of mass and differential coordinates x_, the potential decouples into two terms: 


^'2 


Utot {x+,x-) = k{x+ - -) 


+ + Um{X-). 


(16) 


When considering equilibrium fluctuations at constant trap-to-trap distance A, a linear approximation 
of Um in Eq- (16) may be used. The energy has a minimum at x+ = X, x- = x^. The value of x^ is 
defined by: 


dx- Utot 


= 0 . 


(17) 


Here a;® is the (mechanical) equilibrium distance between the centers of the beads. The position of the 
minimum and the derivative of the potential at the minimum depend on the trap-to-trap distance, so 
that x® is a function of A. In this approximation the equations of motion in Eq. (5) become: 


1+X+ 

y-X- 



(18) 

(19) 


In particular 7 + is the friction affecting the center of mass of the dumbbell, while 7 _ affects the dynamics 
of the differential coordinate. After Bachelor [3] we set: 

7 ;! = (7-i + r-i)/ 2 , 7Zi = 2(7-i-r-i). ( 20 ) 


In brief, 7 is the hydrodynamic friction coefhcient for a single bead, while T is the intensity of hydro- 
dynamic interactions between the two beads. It is important to bear in mind that in principle 7 and 
r depend on the differential coordinate x_. This dependence is negligible in pulling experiments which 
induce a small change in x_, but it becomes clear and measurable if x- is changed to a larger extent, 
e.g. using molecules of different contour length. The equilibrium probabilities generated by (18),(19) are 
given by Boltzmann’s distribution: 

Qeq{x+) = ^exp ^-/3fc(x+- ^)2^ , (21) 

= ^ - {x--x'L) j, 


Peq{X-) 


( 22 ) 
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with normalization factors. The variance of equilibrium fluctuations in a;+ and X- is connected 

to the elastic properties of traps and tether by: 

^2 _ ^2 _ ‘^KbT 

2k ' k + 2km ’ 

Information about the hydrodynamic interactions can be obtained from the time-dependent correlation 
functions of x+ and X-: 


C+{t)^ ((a;+(t)-A)(^^( 0 )_|)) (24) 

C-{t)= ((a;_(t) — )(x_(0) — X® )), (25) 


which characterizes the decay of fluctuations and allows to distinguish the presence of different contri¬ 
butions to the total variance. The computation of the correlation functions yields: 


pc+{t) 

pc-it) 


2k ’ 

2g-(fc+2fcm)7l^< 
k -\- 2km 


(26) 

(27) 


where /? = (KbT) Prom Eqs. (24),(25) the correlation function at time 0 equals the variances given 
in Eq. (23): 

C+( 0 ) = 4 , C_( 0 ) = 4 . (28) 

The hydrodynamic parameters can be retrieved from the correlation functions, as the stiffnesses, k, km, 
are known from Eq. (23) : 


1 1 1 d , , „ , 

t + r = 

1_1 ^ 2 d 

7 r k + 2km dt 

where we used Equations (20),(26) and (27). As discussed in the main text, 7 + is needed to correct the 
error committed by using the force measured in the trap at rest in the Jarzynski Equality (Eq. 17 and 
18 in the main text). 


(29) 

(30) 


S5 A CLOSER LOOK TO HYDRODYNAMIC INTERACTIONS 


In the previous section we have shown how the analysis of thermal fluctuations in a dual trap setup 
allows the measurement of the two scalar parameters entering in the hydrodynamics: 7 “^ and P”^. The 
sum of the two parameters can be obtained through Eq. (29) from the analysis of the time correlation 
function for x^, Eq. (26) . The difference is instead obtained by a similar analysis on the decay rate of 
C-{t). Up to now we have neglected the dependence of the hydrodynamic parameters on the relative 
distance between the two beads. This is justified provided this relative distance does not change too 
much in the pulling experiments. Conversely changing the contour length of the molecule by an order 
of magnitude leads to a larger and detectable change in the hydrodynamic parameters. This offers the 
possiblity of testing our measurements against the theoretical prediction obtained using Stokes equation 
and expressed as a power series in the reduced distance p = {x-)/rb with rj, the radius of the beads and 
(x_) the mean of the differential coordinate. (Note that (x_) > 2rt and thus p > 2). The two quantities 
can be computed to arbitrary precision, taking the first two terms of the series expansion reported in [4]: 


7 


-1 


p-i 


4 (' 4 +°''’"') 
4 (4 4 ''< 4 ^ 


(31) 

(32) 
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We measured fluctuations using ds-DNA tethers of 4 different contour lengths: 8 /rm (24 kbp), 1 /rm 
(3 kbp), 300 nm (1.2 kbp) and 20 nm (58 bp). For each tether fluctuations were measured at different 
forces. The different length of the tethers used in the experiments is such that p assumes values from 
above 6 down to 2.03, near the minimum value 2 taken when the beads are in contact: hydrodynamic 
interactions are monitored from the far-held regime {p ^ 2) to the lubrication limit {p cs 2). The 
experimental data for 7 and F can be confronted with predictions Eq. (31),(32) without free parameters 
as the buffer viscosity and the bead radius are known. 

This comparison is shown in Figure S4, the upper plot compares the measured values for 7 “^ and 
with Eq. (31),(32) (solid lines). For every tether we assumed a single value of p and averaged values 
obtained at different forces. On the force range relevant for pulling experiments, changes in 7 “^ and 
where below 10%. Theory and data do agree, within the error bars for the longer tethers. The data 
for the shortest tether show a deviation from expression (31) and probably more terms of the expansion 
would be needed. Nevertheless in this case the parameters are very close to the contact value, p = 2, 
7 “^( 2 )/ 7 “^((X)) = r“^( 2 )/ 7 “^(oo) = 0.775. The precision of these measurements is limited, in our case, 
by the error with which the radius of the bead is known. One can reduce the dependence of the measured 
value on the size of the bead by using the ratio between the two parameters, which does not depend on 
rt) if not through the definition of p. The results are shown in the lower plot where the measured data 
are compared to those obtained from untethered beads. This allows us to conclude that, at least in our 
conditions and within our experimental resolution, the friction and/or hydrodynamic effect due to the 
presence of a polymer between the beads is not distinguishable. In particular, knowing the value of 7 + 
allows an estimation of the systematic error on unidirectional free energy estimates from single molecule 
pulling experiments based on wrong work dehnitions, as detailed in the main text. 

S6 FREE ENERGY INFERENCE IN ASYMMETRIC SETUPS 

In the case of asymmetric setups P{W) and P{W') do not have such a simple relationship as in the 
symmetric case. In this case a successful inference cannot be based on the FR alone: some additional, 
system-specific, information must be provided. We shall focus on asymmetric systems in the Gaussian 
approximation and use pulling experiments on dsDNA in an asymmetric setup to test our predictions. 
Is it still possible to infer the full dissipation from partial work measurements? The answer is positive if 
we are given some equilibrium information on the system. In this case it is enough to know the trap and 
molecular stiffnesses ks, kA, km i-e. equilibrium properties of the system. But why does direct inference 
fail in the asymmetric case? When discussing symmetric systems we have shown that P{W) and P{W') 
are related by a simple shift: 

P(IU) = P'(IU-A). (33) 

Only the mean of the probability distribution had to be changed, as the variance of the two distributions 
is the same. In that case imposing the validity of the CFR for P{W) yields the unique value of A to be 
used in the reconstruction. This is not true anymore in asymmetric systems. Here both the mean and 
the variance of the work distribution must be changed, which can be achieved by convolution: 

P(1U)a,s = P'*A?(A,S), (34) 

where * denotes the convolution operator and A/"(A, S) is a normal distribution with mean A and 
standard deviation S. Starting from any distribution P'{W') there are infinitely many choices of A and 
S which yield a P(IU)a,s satisfying the fluctuations theorem. Indeed, let us suppose the pair A*, E* is 
such that: 

PA 7 S*(f^) = ^’'*A7(A*,S*) (35) 

satisfies the CFR. Then it is easy to check that Y^(iF2+207/3) satisfy the FT for any (f) {f3 = 

1/KbT as always). In Figure S5 we show the effect of the convolution of P'{W') with different Gaussian 
distributions. In the rightmost column of that plot we highlight three different pairs for which the 
convolved distribution satisfies the CFR. In this situation the inference cannot rest on the CFR alone, and 
additional information about the system is needed. In the general scheme of a Gaussian approximation. 



inference of P{W) is still possible. In the coming subsections we will show that P{W) and P'{W') are 
related by an Asymmetry Factor (AF) given by: 


A.F{hj\^ hm-) — 


~ _ 1 - ks) 

{W)-{W') ~ 15 kA{kB + 2k^y 


where /3 = 1/KbT. The AF is an equilibrium quantity: it only depends on the stiffnesses of traps and 
tether. Knowing the AF allows us to select the unique pair (A, S) such that AF — S^/A and that 
Pa,v,(W) satisfies the fluctuation symmetry. These (A, S) allow the reconstruction of P{W) and thus to 
measure free energy differences or even dynamical quantities like 7 +. The behavior of the AF factor as 
a function of x = kB/kA and y = km/kA is shown in the left panel of Figure S 6 . 

S6.1 PULLING EXPERIMENTS IN LINEAR SYSTEMS 

Free energy inference in asymmetric setups is not as straightforward as it is in symmetric setups. In order 
to recover the full dissipation from partial work measurements we must complement these non-equilibrium 
measurements with some equilibrium information about the system. In this section we will consider a 
Gaussian approximation for asymmetric setups that, as shown in Figure 6 in the main text, agrees pretty 
well with the experimental results. The Gaussian case will be modelled by linear asymmetric systems, a 
class of statistical models for which inference is possible without symmetry restrictions. This amounts 
to choosing Um{x) = \kmX^ in Eq. 1 so that the total potential reads: 

U{yA,yB,X) = ^{VB - VaT + ^Vb + ^ (A ■ (37) 

In order to simplify the following discussion we will now switch to a vector notation. From now on we 
shall denote with X^(t) the vector containing the positions of the traps (remind that B is the reference 
trap at rest with respect to water while trap A moves): 

X^(t) = (At,0), (38) 

and with y the vector of bead positions: y = {yA,yB)- In this vector notation the potential can be 
written in its normal form as: 

^7(y,X^(t)) = i(y-yO) . k(y-y°) -f ^k^gX^, (39) 


where k is the stiffness tensor 


kA km km 
km ks kf] 


fceff is the effective stiffness of the dumbell as a whole: 


^eff 


kAkskm 

det(k) ’ 


and is the vector whose component are the equilibrium positions of the beads: 


y = A-— ^(A),—^(A) 


V det(k) det(k) / 

During a pulling experiment the trap to trap distance will be changed at a constant speed, v: 

Xt = Xo+ vt. 


As a consequence the equilibrium positions of the beads will change in time: 

y"(f) = (a„ + „( - ^(A. + „(), h^(A„ + „()) . 
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During a pulling experiment the free energy of the system changes in time. FRs can be used to determine 
the free energy change in the system from the distribution of the work performed on the system. As 
discussed in the main text, the correct definition of the work in such an experiment depends on the way 
in which the traps are moved. Work is identified as the time derivative of the total potential U: 

W = -V [ kA {yA{s) - aJ(s)) ds (45) 

Jo 

Any different definition of work does not yield reliable free energy differences. An alternative work 
definition can be based on the force measured in the trap at rest (B): 

W' = vf Ub {yB{s) - Xg{s)) ds, (46) 

Jo 

a choice made by many experimentalists. Is it possible to infer the distribution of W from the distribution 
of W' and thus to reliably measure free energy differences? The answer to this is worked out in detail 
below. 


S6.2 W AND W’ AS GAUSSIAN RANDOM VARIABLES 


In order to answer the questions posed in the previous section we need to study the distributions of W and 
W'. In this linear model yA and yB are Gaussian random variables, so that W and W, being linear in yA 
and yB, are themselves Gaussian random variables too. The distribution of a Gaussian random variable 
is determined by its mean and variance. The computation of the mean and variance of the work requires 
the solution of a system of linear differential equations that can be obtained form the Markov Generator 
for the joint process (y, W) or equivalently from the corresponding Fokker-Planck equation. Here by y 
we denote the usual vector containing the positions of the beads while W = (n fA{s)ds, v fB{s)ds). 
We use the vector W from which the random variables W, W' are obtained as; 

W = Wa, if' = -Wb. (47) 

Just as we did when discussing the measurement of hydrodynamic parameters we will use Langevin 
equations to describe the dynamics of our system: 

y =-fik(y - y°) + V, (48) 


where y, is the mobility tensor, describing both friction and hydrodynamic interactions: 



(49) 


and r/ is a white noise compatible with the fluctuation-dissipation theorem, i.e. {rjtVs} = 2/3 ^y6{t — s). 
In this setting the Markov generator L for the process (y, W) is given by; 


Lf(y,W) = 

= (-Mk(y - y°) • Vy + /3-V • - vkD(y - X^) • V^) /(y, W), 


where 


k_D = 


kA 0 
0 kB 


We recall that the Markov generator is the infinitesimal evolution operator in the sense that: 

|(/(y(t),W(f))), = {Lf(y(t),W(t)))y, 


(50) 

(51) 

(52) 


where {)y denotes average conditioned to initial condition (y, 0). Applying the Markov Generator to 
y, (y - {y)y) ® (y - (y).), w, (w - (w)^) ® (w - (w),) ^ and then taking the average, we obtain the 


® j denotes the tensor J such that Jk = (i • k)j. 
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equations of motion for (y)^, cr^ = ^(y _ 0 (y - (y)^;))^, {W)y, C = ((W - (W)y) ® (y - {y)y))y 

and = ((W - (W),) 0 (W - {W)y))y. 


= 

-f^ki{y)y-y°{t)) 

(53) 

d 9 

Id” = 

-2Sym (^pkal + ^ 

(54) 

)(W), = 

1 

1 

(55) 

dt 

-pkC - o-yk^i 

(56) 

d 9 

— 2k£)C, 

(57) 


where Sym(-) in (54) is the operator which gives the symmetric part of its argument. These equations 
are explicitely derived in subsection S6.5 and solved in subsection S6.6, while we shall now comment on 
the result. 

The calculations reported in the next sections show that, after neglecting transients: 


{W{t)) = —v f k£) (y°(s) — X^(s)) ds + k£)k ^k 

Jo 


and 


-2 2 . 
(Tw = 2- - - V t 


The first term on the right hand side of (58) is the integral: 

v^D - (y°(s) - X^(s)) ds^ , 


(58) 

(59) 

(60) 


the difference k£)(y**(s) — X^(s)) gives the equilibrium force in each trap at time s, this is just the 
force that would be obtained at the corresponding trap-to-trap distance A if the pulling were carried out 
reversibly. The two components of the vector kjj ^ (y^(s) — X^(s)) ds'j are of course equal in size but 


different in sign: 


vkn - (y°(s) - X^(s)) ds^ ^ / /i'ev(s)ds,-u J frev{s)ds 


(61) 


The above expression is then just the reversible work: 


Wrev = keS ( XoVt + -V^t^ 


(62) 


where we used Eqs.(39) and (43). The second term in (58) is the dominant non-equilibrium contribution 
(by the way, this is the only term which is linear in time), which arises from the finite pulling speed of 
the protocol. A key element in this dissipation term is: 

Q = kokr^ (63) 


To give some insight on the physical meaning of this quantity we note that p transforms the trap position 
vector X^ into the equilibrium position vector y**: 

yO = pX^. (64) 
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This dissipation term stems from friction the relevant velocity being that of the equilibrium 

positions of the beads. The same element Q appears in the variance. Now, coming to W and W' we 
have: 

(IT) = Wrev + 

9 2 = o 

~ ^^aav t 

iW') = ITrev - ^ABvh 
2 _ 

^W' = 

Based on these expressions we introduce the asymmetry factor: 

2 2 
Ap ^ ~ 

(IT) - {W) ■ 

This quantity relates the distribution of W to W' based on equilibrium information. The AF does not 
depend either on the pulling speed, which is evident from Eqs. (65),(66),(67),(68) or on the hydrodynamic 
parameters. An explicit computation gives: 


(65) 

( 66 ) 

(67) 

( 68 ) 


(69) 


1 AkmjkA - ks) 

P kA{kB + ‘2km)' 


(70) 


S6.3 THE AF CAN BE MEASURED FROM EQUILIBRIUM FORCE TRACES 


The AF can be directly measured from the equilibrium force traces. The equilibrium distribution for 
bead positions follows the Boltzmann distribution with respect to the potential in Eq.(39), i.e. 


P(y) = Z-i exp (-f (y - y°) • k(y - y°)) , 

According to such distribution variance of y is: 

var(y) = 

Moreover, since forces and bead positions are linearly related, f = (y — X^) we have: 

o /o'. _r r~ir _ ^ f k^^kB T km) k^kBkm \ 

pvar^IJ - KCK Ki, - y kAkBkm k%{kA + km) ; ’ 

where /3 = {KbT)~^. Using this formula we get: 

kA = ^var(f)AA + ^var(f)AB 
kB = /3vaT{i)BB + /3vaT{i)AB 
k - var(f)Aj? 

“ *^*''/3det(var(f))- 


(71) 


(72) 


(73) 


(74) 

(75) 

(76) 


These formulas assume that the dynamics of the dumbbell is strictly one-dimensional and takes place in 
the optical plane, i.e. that plane perpendicular to the optical axis. If this is not the case this simplified 
treatment is not valid anymore and out-of-plane fluctuations must be taken into account. A discussion 
of these effects can be found in [1]. Once kA,kB and km are known the AF can be easily computed by 
Eq. (70). The present method for measuring rigidities from equilibrium force traces was applied and 
discussed in detail, in a totally different context, in Reference [2]. We used it in this study to extract the 
values of kA, kB, km for the asymmetric setup. 



12 


S6.4 DERIVATION OF THE EXPRESSION FOR THE AF 

The AF, as given in Eq. (70), can be obtained once the elements of fi are known. We start from the 
dehnition of Q: 

Q = P^/V^P, (77) 


where 


- ^ c-ij, ^ 1 kAiks + km) kskm \ 

^ ^ ^ 'v kAkm ksikA + km) ; ’ 


with 8 = kAks + kAkm + kskm- Here it is useful to switch to a representation in which ^ ^ is diagonal. 
This can be acheived with a 7r/4 rotation: 




Using R we can write: 


V 1 1 


Q = p^RR^/l-iRR^p, 


where, as in Section SI, R^/r ^R is diagonal. The nontrivial contribution which must still be computed 


-/ _ pT- _ ^ f kA{kB + “^km) ksikA + ‘2‘km) ^ 1 ^ B 

^ ^ x/2£: V -kAkm kskA )~ ^8 \-C C 


We can thus conclude that: 


1 (A -C \ ( 7+A 1+B\_ 1 / 7+A2+^_C2 7+AH-7_C- 


28"^ \B C \-l-C j-C J 28^ \ j+AB - 7+H"+7_0- 


By dehnition we have: 


AF = 


2 Qaa — _ 2 A — B _ 1 AkmikA — ks) 

P^AA + ^AB P ^ PkA(kB + 2km)‘ 


S6.5 DERIVATION OF THE EQUATIONS 

In order to derive the equations of motion (53)-(57) we must apply the Markov Generator to the proper 
quantities, recalling that 

|(/(y(t),W(t))), = (L/(y(t),W(t))),. (84) 

We start by deriving Eq.(53). In this case we have that only the hrst term in the generator Eq.(50) 
contributes as the second and third terms involve higher derivatives or derivatives with respect to W. 
We conclude that: 

= {Lyit))y = -pki{yit))y - y°(t)). (85) 

In the same way we can derive the equation for ay, it is useful to recall that 

((y - {y)y) ® (y - (y>!/))y = (y ® y)y - {y)y ® {y)y, ( 86 ) 

so that 

= iLiy®y))y-{y)y ® {y)y ( 87 ) 

we will work out the two terms of the right hand side of (87) separately, starting from the hrst: 

{L (y 0 y))j^ = ((-/tk(y - y°) • Vy + P~^p ■ VyVy) y ® y)y (88) 

The third term in the generator was neglected because it involves a derivative with respect to W. Acting 
with the derivatives on y 0 y)y gives: 

{L (y ® y))j; = -(^k ((y - y°) ® y) - (y <8) (y - yo)) ipk)^)y + 2/3“^A 
= -2Sym ((^k ((y - yO) ® y))^) + 2p-^F 


( 89 ) 
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The second term in (87) is just a time derivative, which can be worked out thanks to (85): 
d 

j^{y)y {y)y = nw?/ - yn » - {{y)y » i Wy - yn) 


;^(y)y o (y)y = -/^k (((y)j^ - yO) O (y)^) - ((y)j, O {{y)y - y®)) = 

= -2Sym {jjk {{{y)y - y°) ® y)) • 


An explicit result fot the time derivative of cr^ is obtained by taking the difference of (89) and (90) 

= -2Sym ((^k ((y - y°) ® y))j^ - ^k (((y)y - y°) ® (y)^)) + 2/?“V = 

= -2Sym {{pk (y ® y))^ - ^k ((y)y ® (y)y)) + 

= —2Sym (yukd^) + 2/3“^/i. 

In the second step above we used the fact that 

(((y - y°) ® y)). = ((y ® y)>, - (y° ® (y).) 

while 

(((y)y - y°) O {y)y) - ((y)^ ® (y)y) - (y° O (y)^). 

Following the order of equations (53)-(57) we should now compute the average of W. This is again a 
simple matter, as the computation proceeds almost exactly as in the case of {y)y, with the only difference 
that now only the last term in Eq.(50) for the generator matters. The result is: 


-(W), = uk^((y),-X^(t)). 


(92) 


The case of C, the cross-correlation of y and W involves some additional computations. The function C 
is linear in both W and y so that only the first and third terms in the generator contribute. The time 
derivative of C = (W — {W)y) ® (y — {y)y) again yields two pieces, as in the case of (87), and we shall 
proceed in a similar way: 

|c = (L(W®y)),-|(IT),®(y),. (93) 

Again we will compute the two parts separately: 

(L(W O y))y = (- (W ® {pk{y - y°))) -k {{vkoiy - X^)) ® y))y = 

= {-pk (W ® (y - yO)) +v{{y- X^) ® y) k^),,, 

and 

/y yy {y)y = -Ak ((W)j^ ® {{y)y - y®)) -k v {{{y)y - y°) O (y)j^) kn- 
Taking the difference of these two terms and using 




(94) 


(95) 


((W ® (y - yO))), = ((W ® y)), - ((W) 


and 


{{W)y o {{y)y - y°)) = ((W)j^ o (y)y) - {{W)y ® yO) 

it is possible to conclude that: 

J^C = -pk{iw - {W)y) O (y - {y)y))y + ((y - {y)y) O (y - {y)y))y^D 
= —pkC -k vdyko- 


(96) 


The only equation left to derive is that for = ((W — {W)y) ® (W — {W)y))y. This computation will 
only involve the last term in the generator: 


|((W - (W),) ® (W - (W),)), = (L (W ® W)), - I ((W), ® (W),). 


(97) 
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As usual we will evaluate the two terms separately and put them toghether as a second step; 

(L (W ® W))y = v{kD (w ® (y - X^)) + ((y - X^) ® W) kD)y, (98) 

I ((W), ® (W),) = vkn (W ® ((y), - X^)) + (((y), - X^) ® W) k^. (99) 

The last equation of motion is thus; 

= +uko((y - (y)y) ® (W - {W)y))y + v{{W - (W)y) ® (y - {y)y))ykD = 

= +2Sym (uk^iC) . 


S6.6 SOLUTION OF THE EQUATIONS 


In the previous section we introduced an average operation {■)y without specifying with respect to which 
measure the average is taken. In this section we will suppose that the initial probability distribution 
for the process is p(y,0) = S{y — y), a deterministic initial condition. The Work is by definition equal 
to 0 at t = 0, again a deterministic initial condion, and the same is true for the different variances and 
covariances beacuse of the initial conditions. The initial conditions for the equations of motion will then 
be; 


(y(o)>!/ = y 

{W{0))y = d2(0) = C(0) = d^(0) = 0. 


( 101 ) 


All the equations we have to solve are linear equations which can either be solved by variation of constants 
or by direct integration. The hrst equation, that for {y)y is a case for variation of constants; 


The solution is 


f (y(o))y = y 

\iiy)y = -yH{y)y-y^{t))- 

{yit))y = (y+ [\ke~^S^is)ds] , 


as it can be directly checked. We recall the definition of y^{t) 

y°(t) = (x + vt- ^^^{Xo + vt),’p^{Xo + vt) 


det(k) det(k) 

In the following it will be useful to write y^{t) — 'tp + where 


, kAkm ksk 

? = r" —^—jwf, V—Twt; 


det(k) det(k) 


In terms of these quantities (103) reads; 


{y{t))y = e / ykef^''‘'{^p + P.t)ds 


( 102 ) 

(103) 

(104) 

(105) 

(106) 

(107) 
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Integrating by parts we get; 


{yit))y = e + ) = 


= I „ -u 






' + e^^‘(V' + Cs)-V--(Mk)“^e^^"d ) = 


- ( y + 


y + e^^‘(V' + ?t)-(Mk)-' - 1 


= I - V-) + (V- + Ci) - (^k)-! (1 - 


= - V-) + y°(t) - (Ak)-' (1 - 


The average (W)j, is given by: 


(W), = k,, f\y{s))y-X'^{s)ds = 

Jo 

= ko ^ -^P)+ yO(s) - (Ak)-i (l - ^ - X^{s)ds 

= ^d( [ (y°(s) - X^(s)ds) - (Ak)“^e“^^*(y - V')! 


- (Ak)-i (s + (Ak)-V^'*)[ej = 

= kzr (y°(s) - X^(s)ds) - (Ak)”' - l) (y - 

+ (Ak)-i ft + (Ak)-'fe-^^‘-l))c) 


The equation for the variance requires some more attention, we first prove a useful identity; let A 
and B be symmetric matrices. The product AB is not in general symmetric 

{ABf = BA. 

As a consequence the exponential of AB, is not symmetric. We want now to prove that: 

( 110 ) 

Using the definition of the matrix exponential we have: 

°° (AB)^ 

A-^e^^ = A“i ^ Y =A-^ + B + BAB + ■■■ = {1 + BA + BABA + ...) A“^ 
fc=o ■ 

= Yi^A-^ = e^^A-\ 


Similarly one can prove: 


(Ae^^) = (Ae^^)^ - e^^A. 


These results will prove useful in soving the equation for ay. We will first solve: 


0-2(0) = 0 


^al = -2^kcr2 + 2|, 
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and then prove that the solution of (113) also solves (54). Equation (113) can again be solved by variation 
of the constants: 


alit) = ^ 


= f 1 _ - 


/? 

-1 


n-1M 


/3 


^ /? 


The original equation for contains Sym (^ka^). Using (110),(112) we can prove that: 


;-i il-i 


/lkd2(t) = Ak (l - _ (i - k^^ = (Akd2(t))^ 


In other words 


^kcry(t) = Sym (^kcry(t)) , 


(114) 


(115) 


so that if cry(t) solves (113) it also solves (54). The equation for C, is again solved by variation of 
constants: 


The solution is: 


C(0) = 0 

4C = -AkC + valkD- 


(116) 


C(t) = ve 




e^'^^ay{s)kDds = 


rt . , _ . 1^-1^ 


= ^ (l - 


D 


/? 




J ^ 

= ue-^^‘(Ak)-i 

= ue-^^‘(Ak)“^ - 2 


/? 


* k ^k£) 
0 /? 
k-iR 


£) 


= u(Ak)“^ (^1 + _ 2e-Mki 

-,-ca 2 k“^k£i 


k-^k 


/3 

D 


a 


= u(Ak)-i (l-e-^^*y 


/3 


As in the case of (IU)y the equation for is solved by direct integration: 


^w = J ‘^v\D^{s)ds = J vkD{^J.k) ^(l-e 


2 Vr^ko 


Id 


= 2v^ / knifik) 1 — 2e + e 


-1 A _ 9„-Mks I -2^iks\ k kp _ 


= 2u^k£)(^k) ^ I s + 2{ij,k) — e 


^ /3 

(Ak)“^ ^_2/iks 


* k-^k,, 

0 fd 


= 2v\d{A)-^ (^t + 2(Ak)-i - l) - - l)^ 

= v^koifik)-^ (2fLkt + 2 - l) - - l)^^ 


k-^k 


D 


(117) 


(118) 


Keeping terms of order 0{t) in Eqs. (108) and (118) we get the expressions anticipated in Eqs. (58) and 
(59) 
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S7 EXPERIMENTS AT LOWER PULLING SPEED 

In the main text we present experiments performed at high pulling speed, in order to enhance the 
dissipation associated with the movement of the center-of-mass. This was done to highlight the fact that 
W does not fulfill the CFR. Nevertheless such violation is still evident at lower pulling speeds. Just 
fpr completeness in figure S7 we show the distributions of IU_, W, W' and IU+ in a pulling experiment 
performed on a dsDNA tether (experimental conditions are identical to those reported in the main paper) 
at a pulling speed of 500 nm/s with A A = 400 nm. 
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CAPTION FOR FIGURE I 

Experimental set-up. Two laser beams, oriented along the z direction, are used to create two optical 
traps. A dumbbell is formed by two optically trapped beads and a molecular tether. The tether is 
oriented along the y direction, perpendicular to the optical axis z. We choose the center of one trap 
(trap B) as the origin of our coordinate system. A denotes the trap-to-trap distance while yA and ys 
denote the positions of the centers of the beads with respect to the reference trap B. 

CAPTION FOR FIGURE 2 

Independent coordinates in linear systems A) Two dimensional histograms of the forces measured 
in a symmetric dual trap setup. The dashed line forms a 7r/4 angle with the coordinate axes and 
corresponds to the definition of /+. B) Covariance (f^ft) as a function of ()»in a symmetric setup. The 
red line denotes 0 covariance (i.e. linear independence). At the bottom of the graph we report the result 
of a 1% significance test, with red meaning dependent and green meaning independent. C) Two 
dimensional histograms of the force as measured in an asymmetric dual trap setup. The dashed line 
forms a 7r/4 angle with the coordinate axes and corresponds to the definition of /+. In this asymmetric 
case /+, /_ do not correspond to the principal axes of the histogram. D) Covariance (f'^ft) as a function 
of (/) in a asymmetric setup. The red line denotes 0 covariance (i.e. linear independence). It is still 
possible to define a coordinate system where the two degrees of freedom are uncoupled, but now the 
definition is different from that given in the main text. At the bottom of the graph we report the result 
of a 1% significance test, with red meaning dependent and green meaning independent. Tests were 
performed on 3 seconds data traces with IkHz acquisition rate. 

CAPTION FOR FIGURE 3 

Independent coordinates in non-linear systems Left panels: force fluctuations of a two-state DNA 
hairpin. Force increases from top to bottom and the hairpin is 1) completely folded, 2) preferentially 
folded, 3) preferentially unfolded, 4) completely unfolded. Right panels: Covariance (/+/f) as a function 
of (j), The red line denotes 0 covariance (i.e. linear independence). At the bottom of the graphs we report 
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the result of a 1 % significance test, with red meaning dependent and green meaning independent. 
Data show that the center of mass (x+) and differential coordinate {x-) are independent. Tests were 
performed on 3 seconds data traces with IkHz acquisition rate. 

CAPTION FOR FIGURE 4 

Direct measurements 7 and F. Left panel: hydrodynamic friction ( 7 , solid symbols) and interaction 
coefficient (F, open symbols) measured from the decay rate of thermal fluctuations. Each point is the 
average upon measurements over 5 different molecules. Solid lines are the theoretical predictions Eqs. 
(31),(32). The horizontal dashed line marks the exact theoretical value at contact {p = 2). Right panel: 
Ratio between the hydrodynamic coefficients (F/ 7 ) as a function of {p — 2)“^. This quantity does not 
depend directly on ri,. Open symbols represent measurements obtained with untethered beads at different 
separations, solid symbols show measurements obtained with tethered beads. The continuous line gives 
the theoretical prediction according to the first two term in the expansion as in Eqs. (31),(32). 

CAPTION FOR FIGURE 5 

Convolution and reconstruction of P{W) from the CFR for the asymmetric case. Each 
row shows the result of a convolution by fixing the value of S and changing A. The top, middle and 
bottom rows correspond to = 0, = 2.5, = 7.5. In each panel we show Pa,e{W) (blue solid 

points), together with Pa,s(— fU) exp(IU) (blue open points) and the experimentally measured P{W) 
(green points). Three different convolutions are found to fulfill the fluctuation symmetry Pa,t.{W) = 
PA,s(~fU) exp(IU), (panels in the rightmost column), but only one of them (bottom right panel) is 
compatible with the AF and matches the true P{W). 

CAPTION FOR FIGURE 6 

Asymmetry factor and inference in asymmetric setups Left panel: The asymmetry factor as a 
function of kA/ks and km/ks- Different curves correspond to different values of km/ks (0.3,0.5,1,2,3 
from bottom to top). The dashed curve corresponds to the limit km = 00 . In the symmetric case 
{kA = = 1) the different curves coincide (S^ = 0). The red point denotes the asymmetric conditions 

in which experiments were performed. Right panel: Inference in the asymmetric case. Different pairs 
(A,S^) yield probability distributions satisfying the CER (blue points and blue line). The asymmetry 
factor (red line) selects a narrow range of possible values, which is compatible with the experimentally 
measured values A = 8.2 iP^T and = 7.2 (A^T)^. 

CAPTION FOR FIGURE 7 

Work measurements at lower pulling speed. The statistics of W, W', W+ and W- are shown. The 
pulling speed v = 500 nm/s is less then half of the lowest pulling speed presented in the main text, but 
the effect on the validity of the CFR for the different work quantities (IF, W', , W~) is still visible. 

500 nm/s lies in the typical range of pulling speed used in single molecule pulling experiments. 
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TABLE SI 


Table S 1: Molecule 1 


Pulling speed 

cov(lT,lT') x^-test^ cov(lT+,lT ) x^Aest^ 

7.2 ^m/s 

4.3 ^m/s 
1.35 /^m/s 

Forward protocol 

0.41 1 0.08 0 

0.35 1 0.09 0 

0.38 1 0.05 0 

7.2 ^m/s 

4.3 ^m/s 
1.35 /im/s 

Reverse protocol 

0.17 0 0.018 0 

-0.05 0 -0.05 0 

0.36 1 0.20 0 

7.2 ^m/s 

4.3 ^m/s 
1.35 /^m/s 

Cyclic protocol 

0.12 1 0.00 0 

0.12 1 0.01 0 

0.34 1 0.09 0 

^ l=Dependent O^Independent 

TABLE S2 

Table S 2: Molecule 2 

Pulling speed 

cov(lT,lT') x^Aest^ cov(lT+,lT ) x^Aest^ 

7.2 ^m/s 

4.3 ^m/s 
1.35 /im/s 

Forward protocol 

0.469377 1 -0.12 1 

0.6170601 1 -0.14 0 

0.5346567 1 0.09 0 

7.2 ^m/s 

4.3 ^m/s 
1.35 /^m/s 

Reverse protocol 

-0.31 1 0.25 0 

-0.37 1 0.20 0 

-0.73 1 0.19 0 

7.2 ^m/s 

4.3 ^m/s 
1.35 /im/s 

Cyclic protocol 

0.26 1 -0.19 1 

0.34 1 -0.18 0 

0.42 1 -0.10 0 


^ l=Dependent O^Independent 



